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ABSTRACT 


The evolution of dislocation sub— structure ahead of a 
moving crack tip with or without the presence of a spherical 
inclusion is examined in mode I cracking using isotropic 
elasticity. This is done by considering the glide force on an 
edge dislocation located at the crack tip moving at speeds in the 
range O < & < c^ and c^< & < c^, where and are the Rayleigh 
wave speed and shear wave speedy respectively. Since the glide 
force on the edge dislocation is determined by the component 
of the crack stress field, the angular dependence of is 
examined at fixed radial vector from the moving crack tip. The 
transformed component of crack stress has maximum value in the 
direction O = 70° for the crack tip velocity less than or equal to 
0.2c . For crack tip velocities in the range 0.2c < & < O.Qc • 
maximum is located in a direction perpendicular to the motion of 
crack tip(d = 90°). For velocities greater than or equal to 
0-8c^, there are two maxima. The first maximum of greater 
magnitude is located along B =90° for all speeds greater than or 
equal to 0.8c^, The second maximum of smaller magnitude is 
located along B = 30° at 0.8c^ and along B = 40° in velocities 
range 0.9c^ < A < c^. For velocities range c^< A < c^ the second 
maximum is located along B = 50°. 

The effect of inclusion on crack-dislocation interaction 
is separately examined. The strain energy density for crack and 
dislocation, and the interaction energy for crack-dislocat ion, 
crack-inclusion and dislocation-inclusion combinations have been 
calculated. The crack propagation characterist ics are analyzed by 
considering the strain energy density of the crack and the 



generation on active slip systems located on the maximum of the 
shear stress lobes is obtained by matching the strain energy 
density of an equivalent bundle of dislocations produced from a 
source to that of the moving crack. This is then incorporated 
into a sub-structure evolution model around a moving crack. 
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CHAPTER 1 
1. INTRODUCTION 

Th« field of fracture mechanics deals with the 
description of a deformable body containing a crack or cracks and 
also characteri zing and measuring the resistance of materials to 
crack growth. Excellent review of the subject had appeared in 
standard text by KnottC133, BroekC483, HellenC393 and Parton and 
MorozorC493 and in a number of ASTM STP Proceedings. The dynamic 
fracture mechanics which is a sub-branch of fracture mechanics 
examines the behavior of moving crack and is applicable to failure 
in high strain rate deformation phenomena like impact loading and 
explosive forming. In these processes crack move at a range of 
subsonic velocities and hence their propagation character ist ics 
are different than at normal velocities. The static equilibrium 
conditions of SriffithC403 and I rwin— 0rowanlI41 3 may no longer be 
valid for dynamic processes. 

There are two ways in which a crack system may become 
dynamic. The first type arises when a crack reaches a point of 
instability in its length where the system acquires kinetic energy 
by virtue of the inertia of the material surrounding the rapidly 
separating crack walls. Such a dynamic state may be realized even 
under fixed loading conditions. The second type of dynamic state 
arises when the applied loading is subjected to a rapid time 
variation , as in impact loading- In both caseSf the deposition 
and or transmission of the energy to the crack tip takes place 
through shock or other mode of wave motion. The stress fields in 
the material is highly concentrated at the crack tip such that the 


externally applied stress is many orders of magnitude lower than 
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that necessary to break the bonds. As the crack propagates in a 
medium which inherently contains a number of microstructural 
inhomogeneities such as inclusions and dislocations^ the following 
interaction arises: (i> crack— dislocation, (ii) crack— inclusion 
and (iii) a combination of <i) and <ii), all of which can alter 
the propagation characteristics of the crack by directly affecting 
the crack driving force. Cracks also act as a important source of 
dislocation generation because of the crack tip plasticity 
associated with a crack. In this work a numerical calculation of 
these interactions between moving cracks, dislocations and 
inclusions is undertaken with possible applications to 
sub-structure evolution and failure in high strain rate 
deformation phenomena. A related dislocation problem and the 
resultant sub-structure evolution under dynamic conditions has 
been examined by ShettyC553. 

1.1 NODES OF CRACK PROPAGATIONt 

In continuum stress analysis for plane cracks three 
basic 'modes' of cracks have been identified. Mode I ( opening 
mode > corresponds to normal separation of the crack faces under 
the action of the tensile stresses; mode II C sliding mode 
corresponds to longitudinal shearing of the crack faces in a 
direction normal to the crack front; mode III Ctecuring mode ^ 
corresponds to lateral shearing parallel to the crack 
front. Extensions in the shear modes II and III bear a certain 
analogy to the glide motions of the edge and the screw 
dislocations. Of the three modes, the mode I is the most 
predominant to crack propagation in highly brittle solids. There 
is always a tendency for a brittle crack to have an orientation 
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that fflinifflizec the the shear loading. Mode I is the most important 
for practical applications, although the other modes and 

combinations of modes I and II do occur. In our discussion, we 
only consider the mode I type crack propagation, 
l.a DUCTILE-BRITTLE FRACTURE* 

The most simple definitions for ductile and brittle 
fracture available in literature are* Fracture accompanied by 
substantial plastic deformation, called ductile fracture and 
fracture without any observable plastic deformation by naked 
eye, called brittle fracture. But at atomistic level plastic 
deformation is also observed in case of brittle fracture. The 
extent of plastic deformation in brittle fracture is smaller as 
compared to ductile fracture. Hence two types of fracture can be 
distinguished on the basis of the extent of plastic deformation 
associated with crack propagation. The more scientific definition 
given by Rice and Thomson Z21 is that a material is intrinsically 
brittle if it can sustain an atomically sharp crack without shear 
breakdown and dislocation emission. A material is intrinsically 
ductile if it is vice versa. A sharp crack is stable in metals and 
other materials that possess disloca'tion sources close enough to 
create a dislocation screen which shields the crack fronts from 
the external stress fields. If the crack is completely 
unstable, then a sharp crack will breakdown under stress by 
continuous dislocation emission and the crack will open with a 
notch of finite angle equal to the angle between a active slip 
pi ane . C 1 ] 

Under dynamic loading conditions the balance between 
ductile-brittle or the cleavage-emission balance may alter. For 
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•xamplcy the ductile— brittle transition in steel depends on the 
strain ratelI133. The recent experiments on ductile copper alloys 
which suggest that cleavage can be generated by stress corrosion 
at the crack tip C123. This is in contrast to the fundamental 
understanding regarding fracture, because copper is completely 
ductile under normal conditions. Sieradzki et <zl.Cll3 have 
proposed that the cracks are generated by the formation of a 
brittle sponge at the crack tip during the corrosion reaction at 
the tip. Lin and Thomson ClOl have proposed that there exists a 
critical velocity above which the driving crack tip stress field 
is sufficient to overwhelm the maximum induced shielding provided 
by any emitted dislocation. Ashby and Embury C93 have suggested 
that in the bcc metals, the response is determined by the nature 
of the force laws at the crack tip. In their view the sufficient 
number of dislocations near the crack tip can be activated to move 
fast enough to intersect and blunt the crack before the crack 
moves on and this leads to an overall ductility. Hart 183 has 
shown that a crack shows a brittle transition at a critical 
velocity, when that crack generates dislocation shielding and the 
inability of the dislocation generation to keep up with the crack 
as the crack speeds up. 

1.3 CRACK-DISLOCATION INTERACTION: 

The brittle fracture in some structural materials is 
followed by plastic deformation and that some crack tip 
deformation occurs during the rapid crack propagation. Question 
arises how the mode of crack tip deformation is controlled and how 
it is related to the onset of brittle fracture. To understand 
these questions, it is necessary to understand the manner in which 
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they interact with cracks during crack propagation. This 
interaction is also expected to determine whether a crack tip is 
to serve as the main source of dislocation. This is important 
because the dislocation emitted from the crack tip shields the 
crack from the applied stress, thereby reducing the local stress 
intensity factor. Dislocations generated from other sources do not 
necessarily shield the crack tip. 

Theoretical analysis of crack-dislocation 
interaction was reported by LouatC71 who used a dislocation model 
for a finite, double-ended, two-dimensional crack interacting with 
a parallel screw dislocation outside the crack in an isotropic 
medium. The same problem for the semi-infinite crack was solved by 
AsaroC351 and Majumdar and BurnsC203 by conformal mapping. Also in 
isotropic medium, Rice and ThomsonC23 formulated the interaction 
between a dislocation of general Burger vector and a crack by a 
thermodynamic technique. In anisotropic media, Atkinson C3] and 
latter Bernett and AsaroC63 presented procedures for solving the 
elastic problem involving a crack and a dislocation. 

The recent renewed interest of dislocat ion— crack 
interaction arises from the realization that a critical stress 
intensity factor may by needed for dislocation emission from the 
crack C23, because of such a critical stress intensity factor, the 
dislocation distribution emitted from the crack in either the mode 
I or the mode II or III situation would be different from those 
obtained previously without such a factorC4-53. The most 
important difference is the existence of a dislocation free 
zone(DFZ) Cl, 143 near the tip of crack. In fact it was the 
observation of the DFZ stimulated the interests in the 
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dislocation— crack intaraction because of the fact that a 
dislocation emitted from a crack shields the crack tip from the 
applied stress, emission process must be intermittent and the time 
and size of the plastic zone depends on the mobility of 
dislocation . 

LouatCTl had made a quantitative eKamination of the 
elastic interaction of cracks and screw dislocation and found out 
that the number of dislocations which can be caused to enter the 
crack and the energy of the crack-dislocation pair. The tendency 
for the crack propagation increases as the number of dislocation 
within it increases. WarrenC361 use plane strain linear elasticity 
theory to develop full field eKpressions for forces on edge and 
screw dislocations in an infinite elastic plane containing an 
elliptical shaped stress free hole and under limiting geometric 
conditions had shown that an edge dislocation with Burger vector 
normal to the crack front, moved away from the crack tip and those 
with Burger vector parallel to the crack front, accumulated at the 
crack tip. 

1.3.1 DISLOCATION-FREX: ZONE CDFZOs 

Bilby et al.tAl found out dislocation distribution 
emitted from the two ends of a double-ended crack in the mode II 
or III situations and shown that there was no barrier for 
dislocation emission. As a result, when the crack had emitted all 
the possible dislocations, the stress intensity factor at the tip 
of the crack becomes zero. Using similar criterion, the emission 
of four arrays of edge dislocation in the mode I was studied by 
VitekClSl using the interaction of an edge dislocation with a 
crack and by RiedclCSl using the dislocation model for a crack. In 
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all these cases the dislocation density approached infinity at the 
crack tip. 

The Bilby, Cottrell and Swinden model (BCS) C43 dose not 
allow crack growth caused by shielding of dislocation pile-ups at 
crack tip. The problems associated with the BCS model are due to 
the absence of an energy barrier for dislocation emission from 
crack tip. The Ohr and ChangC16yl73 assume a DFZ near crack tip 
and found out a new distribution different from the BCSCC41. 
Between the crack tip and plastic zone, there is a region free 
from dislocation, which is referred to as the dislocation free 
zone. In this new distribution, the distribution function started 
from zero at one end of the DFZ increases to a maKimum and 
decreases to zero again at the end of the plastic zone. The stress 
intensity factor at the crack tip was no longer zero and depends 
on the length of the DFZ The Ki at the crack tip would be zero 
only if the DFZ was absent. The same type of solution was found by 
Majumdar and Burns C183 for a semi-infinite crack. Rice and 
ThomsonC2!l in analyzing the dislocation-crack interaction 
suggested that there should be a barrier to dislocation emission. 
With a barrier the equilibrium distribution is such that there is 
a DFZ at the crack tip. One such barrier is the image stress on a 
dislocation. 

The DFZ near crack tip was observed byCl,14] in 
transmission electron microscopy (TEIi) . Chia and BurnsC341 have 
observed distribution of dislocations ahead of a crack tip in bulk 
single crystal of LiF using the etch pits technique under mode I 
loading condition and found that the maximum density of 
. dislocations was located approximately 6-10 //m in front of the 
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crack tip. This implies that there is a DFZ in the immediate 
vicinity of the crack tip. This observation is in agreement with 
the results of in situ TEM studiesCll. 

1.3.8 DISLOCATION EMISSION! 

The direct observation of dislocation emissions from 
crack tip are confirmed by Kobayshi and OhrCll, who did in 
shear cracking (mode III) experiments on their foil of Mo and W, 
Ohr and NarayanC191 on similar experiments in stainless steel and 
Majumdar and Burns C203 who showed dislocations at crack tip in 
bulk MgO revealed by TEM after thinning. The dislocation 
distribution near the crack tip revealed a region free of 
dislocation next to the crack and this was attributed to The 
attraction between the crack and its emitted dislocat ionsC203 . 

Chia and BurnsC343 in their experimental etch pit 
studies, under mode I loading conditions, choose the crystal 
orientation such that the cracks were made to propagate along the 
<100> as well as the <1101 planes. Experiments reveal that the 
edge dislocations were emitted on two symmetrically inclined slip 
planes. The slip planes intersect the crack plane along the crack 
front. The Burger vectors of these dislocations had a component 
normal to the crack plane and hence these dislocation blunted the 
crack tip during emission. The second type of crack tip geometry 
observed consists of arrays of screw dislocations on two mutually 
perpendicular slip planes that intersected the crack plane at 45? 
They concluded that crack tip was again the source of the 
dislocations. These dislocations intersected ahead of the crack 
tip and formed sessile dislocations of a <100> type, which could 
impede crack propagation. These dislocations do not blunt the 
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crack tip but -rather cause the crack front to jog. 

Rice and ThomsonC23 proposed that if the applied stress 
is sufficiently high so that the attractive interaction is less 
than or equal to the core radius (ro) of a dislocation then the 
emerging dislocation will always be on the repulsive side of the 
force distance curve. The dislocation is emitted spontaneously in 
this situation. In addition, as suggested by Bilby et at.C43, the 
force on the dislocation must be greater than the lattice friction 
force, bo"^, in order for the dislocation to glide away from the 
crack. The dislocation emission condition proposed by ThomsonC21II 
and Lin and ThomsonC22] is that the combined stress intensity 
factor (included all the three modes), should be greater than or 
equal to some critical stress intensity factor for dislocation 
emission. For a moving crack number of dislocations emitted at 
the steady state is smaller than the equilibrium number emitted 
from a stationary crack. The faster the crack moves, the smaller 
is the steady state number. Also for a moving crack there are 
saturation numbers even at zero lattice friction (no saturation 
for stationary crack) and with increasing lattice friction, the 
effect of crack motion diminishes. The moving crack also has a 
smaller plastic zone size (the distance between the farthest 
dislocation and the crack tip ). After each emission, the crack 
stops momentarily because of a sudden reduction of Ki. After the 
emitted dislocation moves away, the crack speeds up again until Kz 
increases to some critical value C571. 

In mode II and III, edge and screw dislocations 
respectively, are emitted on slip planes that are nearly coplanar 
to the crack plane. In these instances, the crack tip remained 
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sharp, since Burger vector is parallel to the crack plane. Cracks 
often propagation in a zig-zag (nannerC193. In this process a crack 
propagating along a path suddenly stopped. It could be seen • that 
many of the dislocations emitted from the crack tip are piled up 
ahead of the crack tip in a con^lex tangle thus making the motion 
of subsequently emits dislocations difficult. The crack remains 
stationary for a while and then emitted dislocation of a 
different Burger vector on an inclined slip plane. After a number 
of dislocations are emitted, the crack turns and propagates along 
this second slip direction C193. 

Zhang and LiC333 clarified the confusion in the 
literature about the interaction between a screw dislocation and a 
finite crack. It makes a difference whether a the dislocation is 
emitted from the crack or comes from else where. For example, both 
the image force and the shielding effect are always larger for a 
dislocation emitted from the crack than for the same dislocation 
originating elsewhere. They become equal if the crack is very long 
(semi-infinite) . 

1.3.3 DISLOCATICW SHIELDING OF A CRACK TIPt 

In addition to dislocation emission at the crack tip, 

dislocations also play an important role in material toughness 

through their strong elastic interaction with cracks. w»hen a 

dislocation is emitted from a crack tip it moves away from the 

crack tip because the net force on it is posi ti veI323 . The force, 

however, decreases as inverse square root of distance from crack 

tip (-■,-) and the dislocation comes to rest where the crack force 
rr 

is balanced by the lattice frictional force. The frictional force 
which is termed as the deformation resistance, is usually a 



fraction (<1) of the flow stress of the materialsCSOl . The region 
between the crack tip and the area where the dislocation comes to 
rest is dislocation free and thus the physical origin of the DFZ. 
The presence of a dislocation in the vicinity of a crack tip 
modifies the stress field of the region because the interaction 
between the crack and the dislocation. This phenomenon is very 
similar to that of work hardening in plastic deformation. During 
plastic deformation, a dislocation generated from a source enerts 
back stress on the source. Due to back stress, the stress at the 
source is reduced accordingly. In order to operate the source and 
hence to continue deformation, the applied stress must be 
raised. This increase in the applied stress is one of the causes of 
work hardening observed during plastic deformation. 

The stress field in the immediate vicinity of the crack 
tip is reduced in the presence of dislocationsC31 3 . this decrease 
in the stress field reduces the local stress intensity factor K, 
which is less than the applied stress intensity factor C301. This 
reduction in the local stress intensity factor caused by the 
presence of a dislocation is referred to as dislocation shielding 
of the crack tip. 

1.3.4 CRACaC TIP BLUNTINGt 

After Kelly et al.C293 it is generally believed that 
whether a solid is ductile or brittle depends upon whether the 
ratio <ymax/Tm<3oc is large or small. Here ermax is its theoretical 
tensile strength and ttoox is its theoretical shear strength. If 
this ratio is indeed a critical material constant that determines 
to a large extent the ductility of a solid, then an atomically 
sharp crack might be expected to be mor-e easily blunted the larger 
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is this ratio. Ulhen the ratio owax/Tinax is equal to or larger than 
about seven a mode I crack tip should blunt before brittle crack 
propagation occurs. The fee metals satisfy this condition. The bcc 
metals Iron and Tungsten have O'nKxx/Tmax ratio smaller than seven 
and thus cracks in these metals not blunt. 

WeertmanC283 used a smeared-out dislocation type 
analysis to determine condition for the blunting of an atomically 
sharp crack that is open in tension. When the ratio of the 
theoretical tensile strength o'max to theoretical shear strength 
rmox is large, pairs of dislocation of opposite sign are created 
on an inclined slip plane whose intersection with the crack plane 
coincides with the crack tip. Dislocation of one sign move into 
the lattice away from the crack tip region and dislocations of 
other sign pile-up against the crack tip Csee 1.1a J>. The 
barrier to the pile-up is the surface energy of the new surface 
that created during blunting. If the stress is raised further, 
enough piled up dislocations are created to enable the surface 
energy barrier to be penetrated. The piled up dislocations enter 
the crack tip and in doing so blunt the tip (see l.lb>. 
Wee-rtman concluded that perfect dislocations are not emitted 
directly from the crack tip as assumed by -Rice and ThompBonC23, 
they assumed that a perfect dislocation is emitted directly from 
the tip itself onto an inclined slip plane. The dislocation move 
away from the tip, as shown in Fig. (1.1c) and blunt the crack. 

The effect of impurity atoms in the crystal lattice, the 
pinning effect of impurity atoms make it more difficult for a 
dislocation to move into the lattice away from a crack tip, the 
qualitative effect of impurity atoms is to make it more difficult 
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■for a crack to be blunted. Impurity atoms in an atmosphere that 
can the crack surface would decrease the surface energy of solid 
and thus decrease the ratio o’max/ rmax . This effect also makes 
brittle propagation of a crack more likely. 

1.4 POINT DEFECT INTERACTIOMj 

The interaction between second phase particles or • point 
defects with crack and dislocations arises due to misfit stress 
field present around the second phase particles. The misfit stress 
result from the second phase particle not exactly fitting into the 
hole in which it sits. For a spherical inclusion Motto and 
NabarroC273 showed that such stresses are proportional to the 
volume of the second phase particle and decay as the inverse cube 
of the distance from the second phase particle. The distortions 
produced around the spherical inclusion are spherically symmetric 
and only normal stress and strain field are present . 

1.4.1 DISLOCATION-POINT DEFECT INTERACTION: 

The hydrostatic stress field associated with an edge 
dislocation can directly interact with an inclusion producing 
spherically symmetric distortionsC44D . Whereas this interaction 
dose not exist in the case of screw dislocation because of the 
non-existence of normal stress components, as they can provides a 
strong interaction if the distortions around inclusion are 
asymmetricCSl 3 . The interaction between a dislocations and a 
second phase particle has three distinct contributions, arising 
from: 

<i) Mismatch in elastic constants between matrix and 
second phase particle. 

(iil Stresses resulting from misfit, and 
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(iii) The concentration cf applied stress by the second 
phase particle. 

If — M — > 1 , the dislocation is repelled from the second phase 
particle and if — ~ — < 1, the interaction is attractive^ where 
and S is the shear modulus of inclusion and matrix, respectively. 
The interaction between moving dislocations and inclusions, may 
result either in the formation of a dislocation loop around the 
inclusion or cutting through the inclusion depending on the 
characteristics of the particle. For hard inclusions, formation of 
dislocation loops are most likely and for soft inclusions, 
dislocation cuts through the inclusion causing slip. 

1.4. a CRACK-POIMT DEFECT IKTERACTION: 

There are a number of ways in which a point defect 
interacts with a cracks a first order size effect arises from 
interaction between the stress field of the crack tip and the 
atomic displacement around the point defect. There will also be a 
second order size effect due to non-linear strains in the vicinity 
of the crack tip and the importance of this interaction will 
depend on the magnitude of the third order elastic constantsC553 . 
Since in elastic theory we represent a point defect by an elastic 
inclusion there will an induced inhomogeneity interaction due to 
the fact that the material in the inclusion responds to crack tip 
field under a different Hooke's law from the surrounding material. 
LiuC263 had discussed the strain field interaction between a 
crack, loaded in tension and a point defect. This interaction is 
of importance since the gradient of the interaction energy 
provides the driving force for the drift flow of point defects to 
the crack tip. EshelbyC373 had shown that the second order 
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inhoroogeneity effect attracts soft inclusion to the crack tip but 
repels the hard inclusions. This is just opposite the to the first 
order size effect. 

The second phase particles located in the near tip field 
of a propagating crack perturb the crack front, causing a 
reduction in the stress intensity factor. The reduced intensity 
depends on the character of the particles and the nature of the 
crack interact ionCS4!l . Two dominant perturbations exists, termed 
crack bowing and crack deflection. It is essential to discriminate 
between crack deflection, which produce a non— planar crack and 
crack bowing, which produces non-linear crack f rontC24,253 . Crack 
bowing results from resistant second phase particles in the path 
of propagating crack, causing the stress intensity Ki along the 
bowed segment of the crack to decrease, hence increase toughness 
of the material. Deflection toughening arises whenever interaction 
between the crack front and the minor phase produces a non— planar 
crack, subject to a stress intensity factor lower than that 
experienced by the corresponding planar crack. The non-planar 
crack arises from elastic modulus mismatch or misfit strain. The 
sign of misfit strains determines the direction of def lectionC233 . 
Specifically a second phase particle with a greater elastic 
modulus than the matrix produces a tangential compressive strain 
near the particle-matrix interface and divert the crack around 
particle. If elastic modulus is lower than the matrix, induces 
tangential tensile strains, causing the crack to deflect towards 


the particle. 
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CHAPTER a 

a.l INTERACnOM BETWEEN A MOVING CRACK AM> DISLOCATICWs 

Suppose that a mode I crack is extending in a 
homogeneous and isotropic elastic body in such a way that the 
state of deformation is two-dimensional plane strain. A 
rectangular (x,y,z) coordinate system is introduced so that the 
crack edge is parallel to the z-axis and the displacement vector 
of each particle is parallel to the x,y-plane, as shown in 
Fig. (2.1). The crack grows in the x-direction with a uniform 
velocity so that the coordinates of the crack tip in the 
x,y-plane are x®a(t), y=0, where 'a' is a continuous function of 
time. The crack tip speed is restricted to the range 0< & < c^ 

and to be continuous. Two types of waves are originated from the 
crack tipC381t 

<i) The dilatational wave with the characteristic wave 


speed c^. 

(ii) The shear wave with the characteristic wave speed c^. 

The dilatational and shear wave speeds for an isotropic elastic 
material are defined in terms of the mass density p and the shear 
modulus G and Poisson's ratio v The relationships are 


C39]! 


1 - V 

1 - 2 vj 


( 2 . 1 . 1 ) 


c^»[— (2.1.2) 
The two types of waves are independent in the interior of a solid. 
At boundaries such as the faces of a crack for example the two 


types of waves interact. 

2.1.1 MOVING CRAC3C TIP STRESS COMPONENTSs 

The X and y axes are fixed in space, so that within the 


small displacement assumption. 


a 


local rectangular coordinate 
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system If and C ^ (£ = 1 or 2) are considered at the crack tip. In 
terms of the position x*a(t) and the speed k of the crack tip, the 
position with respect to local coordinate system is given by the 
following relations i 

K ~ X - a(t) 


and , 


Where , 




i = l,2 


*■■[■ - (-fjf 

V [■ - (J;)'] 


t/z 


(2. 1.3> 


A local polar coordinate system is also introduced, where 


and , 


y( ? cp 




<2. 1 .4) 


Using relation (2.i.3>, the equation <2.1.4) reduce to the 
following form : 


r^= r ^1 - (1- <5*) sin 


2 ^ 1 1 /^ 


tand . » <5 . tan© 

i i 


<2. 1 .5) 

For a detailed solution see appendix (A). 

For a state of steady growth, the crack tip stress field is given 
by the following equations C393 : 

Ki 




f , . (©,4) 


( 2 . 1 . 6 ) 


■/(*nr ) ij' 

Where Ki is stress intensity factor. The function f^^,(©,i> 
represents the angular variation of stress components for any 
value of crack tip speed k and are given below : 
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fii = .5^) - 4.5 a <=os(e4/4> 


yr. 


r V 7- 

12 r 


f22 = 


fl2 * 


_ ^ f r< . jt* 'k* cos(di/*> , cos(82/2) 

"Er[U'^ ‘^2 J — ^ 




2 ( 1 + r «in(di/ 2 ) sin ( 82 / 2 ) 1 

[ ^ y j 


fas » V (fli + f22) 


(2.1.7) 


VIhortt , 

D » 4 <5^6^- (1 + <5*)® ( 2 . 1 . 8 ) 

[1 - (1- <5*) sin*8 ]*"'* (2.1.9) 

The components of stress have an inverse square root dependence on 
radial distance from crack tip and a characteristic variation with 
angular position around the of the crack for crack tip speed in 
the range o < & < c^. In equation (2,ls8)^ D o as where 

c^ is the Rayleigh wave speed. Due to the appearance of the 
factor D in the expressions for these fields^ however^ the 
algebraic sign of the square root singularity depends on whether 
the crack tip speed is less than or greater than the Rayleigh wave 
speed of the material c^. The crack faces separate behind the 
advancing crack tip and eqns(2.1.7) shows that the stress 
components have two-fold symmetry and the traction on the 
prospective fracture plane ahead of the tip is tensile if o < 3 l < 
c but that it is compressive if c^< k < c . This feature implies 
that there is a net flux of energy out of the advancing crack tip 

for c < & < c 1 where as there is a net flux of energy into the 

a z 

crack tip for crack tip speed less than the Rayleigh wave speed c 



C403 


2*1*2 STATI<»<ARY CRACK TIP STRESS COMPONENTS* 

It can be verified that the equations{2. 1 .7) do indeed 
reduce to the equivalent results for a stationary crack in the 
limit A->o. I rufin— Wi 1 1 i ams stress fields for a stationary 
crackC41 3 t 


Where , 



Ki 

’TTznr) 


f . Ae) 

tj 


( 2 . 1 . 10 ) 


fti(^) » cos(d/*)£l - sin (d/*>Bin <3d/z)3 
fzz(0) * ca5(d/z>[l + sin (d/*)sin (3^/2) J 
f33(d) * V [fli<d) + f22(d)] 

ft2(d> * cos (^/2)sin <d/2)cos ( 3 ^/ 2 ) (2.1.11) 
Some care however is required in arriving at the equations (2. 1. 11 ) 
from eqns(2.1.7> as & o (see appendix B) . 


2*1*3 STRESS FIELDS OUTSIDE A SPHERICAL INCLUSION* 

Common Inclusions such as A1 0 . tInO and 5i0 are 

2 3 ' Z 

present in many structural materials as they get introduced during 
manufacturing. These inclusions act as sites for nucleation of 
defects like cracks and voids and affect the ductile brittle 
fracture characteristics. Inclusions located in the near tip field 
of a propagating crack perturb the crack front, causing a 
reduction in the stress intensity factor. Two dominant 
perturbations exist, crack bowing and crack deflection as reported 
by C423. The stress components are obtained by assuming that a 
spherical inclusion of radius r» occupies an interstitial site of 
radius r' in an original lattice and distortions produced around 
the inclusion are spherically symmetric and only normal stress and 
strain components are presents. Eshelby C433 had given the 
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outside stress fields around an isotropic inclusion in an 
isotropic medium as: 


O' * -2o’ 
rr yf \ff 



®_ii 

3 


f \ V f r* ^3 

( t - t. ) (— F— ) 



The crack-dislocation geometry considered to examine 
force on the dislocation is shown in Fig. <2. 2). The slip plane 
intersects the crack plane at an angle 9 along the crack front. The 
dislocation is parallel to the crack front and thus the geometry 
is appropriate for treating a dislocation that is emitted from or 
absorbed at the crack tip. The force on the dislocation is 
composed of three terms: 

1) The force on the dislocation due to the crack 


stresses. 

2) The image force on the dislocation from the crack 
surfaces, and 

3) The interaction between these dislocations and other 
dislocations present in the matrix. 

The image force and the interaction between the dislocation and 
other dislocations are negligible in comparison to the dominant 


crack-dislocation interaction. 


The image force which is always 
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attractive, however, responsible for the formation of the 
dislocation-free zoneCDFZ) near the crack tip and is inversely 
proportional to dislocation distance from the crack tip 1321. 

The transformed crack stresses, a‘ , acting on the slip 
plane of dislocation obtained by using the following rotation 
matrix : 

coe& sin& 0 
-sin^ cos$ 0 

001 (2.1.14) 
The components of stress are given byt 


^j'l 


(2.1.15) 


Mhere, is the crack tip stress components and the direction 


cosines a are given by (2.1.14) ast 


a » a a cos6 
11 22 

01 * “ «_» sin© o « 1 
tt Zi ss 


\ I.T.. 




tK.xi2m 


<n* <x » a a a . o 
9JL 19 29 92 

The resolved shear stress <y' acting on the slip plane, 

1 2 


(2.1 . 16) 


' "i/t Si Si 


<x fa cr+ao’+ao’l+a rao'+ao'+ao*! 
llL'**i'^ll ^ZZ 12 28 laJ 12 21 21 22 22 28 aS-* 

+ a Ta<r+acr+a<yl 

IsL 21 ai 22 92 23 83-> 


Using (2.1.16), <7' becomes 

2 


cr' a C09Z& a + -ST" *in2$ O' ^ ^ 

12 12 2 22 11 


(2.1.17) 


For stationary crack, using (2.1.9) and (2.1.10), b9com«s 


a y sin© cos (d/2) 

12 2r2tTr 


(2.1.18) 


For moving crack, using (2.1.6) and (2.1.7), o'' becomes 




xzonz& f + -=r— sin2^ 

12 ^ 2 2 




(2.1.19) 


The glide force on an edge dislocation that lie on a 
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slip plane, slip plane intersect the moving crack front, is given 
by following relation: 

F = cy'^ b (2.1,20> 

Where F is the force per unit length of dislocation and b is the 
Burger vector of dislocation . 

a. a EFFECT OF INCLUSION ON CRACK-DIS-OCATION INTERACTK»lt 

First, we separately examine the interaction between the 
crack and dislocation and dislocation and inclusion and then by 
taking all the three crack, dislocation and inclusion together. . 
a,a.l CRACK-INCLVSION INTERACTION: 


As a first approximation, we treat inclusion as a 
spherical particle that fits into an undersized or oversized hole, 
resulting in spherically symmetric normal stresses. Such type of 
particles interact strongly with the hydrostatic stress at the 
crack tip (see Fig. 2. 3). 

The hydrostatic stress at the crack tip is given by the relation: 

- —if • 

Following Cottrell's calculationC443 one obtains the elastic 
interaction: 

- —if <7. . . AV (2.2.1) 

Where AV « An£ r'* is the volume change of the spherical hole. 

For moving crack, using (2.1.6) and (2.1.7), the elastic 
interaction energy Urn given by following relations 

2 N /’r a 




u» j fi * ccce./. 

3D rtrzr Y r ^ z ^ ^ 


( 2 . 2 . 2 ) 


Urn 




t (2.2.3) 

Similarly, for stati on ary c r ack , us ing(2.1.^) and (2.1.10), the 


elastic interaction energy Us are : 
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Us = - 


Sn.Ki. (1+ * 


3y'2nr 


cos ( d/ 2) 


(2.2.4) 


M Ki 


cos (d/ 2 ) 


(2.2.5) 


Where y M * 


Qrr. (1+ v) .s.r‘^ 

3 Vzn 


is e constant. 


The equipotential curve for moving and stationary crack have been 
drawn in Fig. (2.3) and Fig. (2.4), using (2.2.3) and (2.2.5) 
respectively. The equipotential curves are the contours having 


constant values of 


“dW <5*)-COS(dl/2) 


for moving crack, for different normalized crack tip speed. For a 

•”i/ 2 

stationary crack, the contours have constant values of r cosd. 

SLi 

The negative of first derivative of interaction energy( - ■■ ) 

gives the force on the inclusion. This force is attractive if £ 
is positive and repulsive if £ is negative. 

2. 2. a DISLOCATION-INCLUSION INTERACTIONx 

For the case of a straight positive edge dislocation in 
elastically isotropic materials, the stresses in terms of the 
three orthogonal coordinate axes, are given by the following 


equations C4S]t 




b Bind (2 cos ) 


d' ■ 


b Sind cos2d 


d' . 

O' » - to 

as 


2 b V sind 


d' . 

O' ® to 
12 


b coad cos2d 


( 2 . 2 . 6 ) 


Where , 


to * 


2n (1 - V) 


The hydrostatic stress of a positive edge dislocation at position 


(r,0) (Fig. (2.5)) from the dislocation is 
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. = - 4 - ) 

m 3 tl 22 99 

(1+ v) 6 b sin6 


(2.2.7) 


3rt (1- v>> r 

Here, « 9 ain the inclusion is treated as a point defect causing 
spherically symmetric distortions. Hence the interaction only 
occurs with the hydrostatic component of the dislocation stress 
field. The elastic interaction energy between the dislocation and 
the inclusion is obtain using (2.2.1): 


Ut 


4 (l-»- y) S b r* sing 

Zn <1- v») r 


Where , 


Ut 


^ sin^ 

A £ 

r 


( 2 , 2 , 8 ) 


^ _ 4 (1+ V) 6 b r'* 

* 3it (1- V) 

hk>te: In section <2.2,1^ and <2.2.£>» the interaction energy 
inclvtdes only the energy external to the inclrxsion. It neglects 
the strain energy change within the inclusion. If the energy 
change of the inclusion is includedr the interaction energy is 
modified by a constant factorial . 


2.2.3 CRACK-DISLOCATION-INa.USION SYSTEHt 


The geometric configuration used in our calculations is 
shown in Fig. (2,6), We assume that the nature of each source of 
ctpcss is unaffected by the presence of the others. So that we can 
superimpose, stresses from different sources over each other. The 
total shear stress experienced by the dislocation is obtained by 
the algebraic summation of stresses from the individual sources, t 
e. the crack and the inclusion. 

The dislocation and the inclusion are taken in the same 
orthogonal coordinate system(x^ ,y* ,z' ) . The shear stress on the 
(dislocation due to the inclusion is obtained by using the 
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following transformation matrix: 


sin^ cosyf 

-siny cosip 

cosy 

sin^ siny^ 

cosy/ cosip 

siny 

cos^ 

0 

— sin^ 


<2.2.9) 

Using •qna(2.1.11) and <2. 1.14). one obtained the for 

inclusion as below: 

^ a a ^ (2.2.10) 

iz 11 zi rr iz xz ^ ^ 13 z» ^ ^ 

Using (2.2.9) y (2.2.10) becomes on simplification: 

^ V ) (— F“)^ ^ (2.2.11) 

Nowy total resolved shear stress on dislocation, obtained by 
using (2.2.11) and (2.1.17), for stationary crack: 


0 ^'* » <?*• + cr* (2.2.12) 

iz iz - IZ 

and, for moving crack, using (2.2.11) and (2.1.18) : 

= <r'"‘ + </ (2.2.13) 

iz iz - IZ 

Here superscripts and 'm' refer to stationary and moving 

crack, respectively. The ± sign put because e may be positive or 
negative. 

The force on dislocation, per unit length of dislocation is given 
by : 

F » cr'*’ . b (2.2.14) 

iz 

-2 

For a typical valu# of r»/r » 0.8, c - +0.55x10 and ju = 156Gpa 
for A1 O in Fa matrix C473 and v - 0.33, Ki * 65fipaym, r » lOiLim 

and G » 736p* for Fe, the force on dislocation have been 

calculated for different normalized crack tip velocity and plotted 


T.ni , m ^ I 

12 12-12 


in Fig. (2.6) • 
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a. 3 STRAIN EVERGY AND INTERACTION ENERGY OF CRACK, DISLOCATION AND 
INCLUSION SYSTEM* 

2.3.1 STRAIN ENERGY GF A CRACK* 

The extension of crack requires the formation of a new 
surface, with its associated surface energy. The elastic strain 
energy density stored in the body is released in the process of 
generation of new surface. The strain energy density is given 


2 


Where , 


1 + V 


<?. . - 


°'hk 


<2. 3.1) 


1 when £=j' 

* 0 when 


(2.3.2) 


For mode I loading, using relation (2.3.2) and (2.3.1): 


* 2 


+ <y* + - 

E t 11 at sa^ 


fo' 0 ’+ 0 ’ 0 '+ 0 ’ 0 'l + ■ or 

E ^ It 22 22 as ss aJ 2G ±2 


1 2 


Now, for a moving crack using relation (2.1.6) and (2.1.7), the 
above relation becomes: 


1 1*" - M 
c * 2nr 


r-^A=-rt*+ f* + f* 1 ^ ff f+ff+ffl 

L 2 E V 11 tz smJ E ^ 11 22 22 as aa 11-^ 


Whcra , 
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U = f(r).f(d,i) 
c 


f (r> 


(2.3.3) 


Here, f(r) is a function of radial distance 'r' and f(G, ) is a 
function of angular displacement G and crack tip velocity 'A'. 
Similarly, for a stationary crack using relation (2.1.9) and 
(2.1.10), strain energy density is given by the relation : 


U = f(r).f{e) 
c 


(2.3.4) 



'll 


Where » 


f (r> » 


Ki‘ 


2irr 


fo) --sV ♦ ’L'®’ * ^»'®0 - -r 

f^<e)t^(0) - f„<0)f„<e>) * 

The enguler variation of strain energy density for moving and 
static crack have been plotted in Fig. <2. 7). 

For the volume element r d® dr dz, the total strain energy per 
unit length of crack front is given by the relation i 

(2.3.5) 

Since there is a singularity at crack tip, the lower limit is 
taken as some finite small value a*. 


X*, U r d© dr 

C O R» C 


For moving crack, using (2.3.5) and (2.3.3), we get 

e” . -- «•’ /«" f(e,t) d« 

Similarly, for stationary crack, using (2.3.4), we get 

E* . de 

c 2 n o 


(2.3. £>) 


(2.3.7) 


The integral (2.3.6) and (2.3.7) have been solved with the help of 
Simpson's three point numerical integration formula. Numerical 
results have been plotted for different crack tip velocity in 
Fig. (2.8) . 


2.3.2 STRAIN ENERGY OF AN STATIC EDGE DISLOCATION: 

The total strain energy of a dislocation is the sum of 
the elastic strain energy plus the energy of the core of the 
dislocation. As a good approximation, the core energy is added to 
the elastic strain energy by taking core radius r^» b, so that the 
total strain energy per unit length of edge dislocation is given 


byC453 i 


U 


S b 


4 n 



(2.3.8) 
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a. 3. 3 STRAIN ENERGY DENSITY OF A >K>VIHG DISLOCATIOMt 

We consider e subsonic discrete dislocation in uniform 
motion in an isotropic media* The moving straight dislocation is 
considered to remain straight during motion referred to its own 
cartesian frame and lies parallel to the a— axis and moves in the 
X'-direction with velocity v within the subsonic velocity range O < 
V < c^* The dislocation source producing the dislocation lies on 
the slip plane that intersects the crack front. Let the source 
produce an equal number of positive and negative edge dislocations 
at any time as a result of the moving crack stresses. The total 
number n of dislocations matching the crack strain energy density 
to that of a giant dislocation of Burger vector nb required to 
replace the moving crack stress fields on the slip plane of 
dislocation. The stress field is then obtained by considering a 
giant dislocation of Burgers vector nb present on the slip plane 
as shown in Fig.(2.15> at all instant of crack motion and then 
breaking down to n elementary dislocations after the passage of 
the crack. The stress field associated with a uniformly moving 
edge dislocation of Burgers vector b is given by the following 
rel at ionsC431 : 


p s 
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P„ = v(p 


+ p 1 
11 > 


(2.3.9) 


Where , 


-[■- firl 


1/2 


1/2 


' ^ t )'] 

The strain energy of a dislocation is proportional to the square 
of the Barbers vector of the dislocation. The total strain energy 
of the giant dislocation of Burger vector nb is going to be 
considerably greater than the sum of the strain energy of n 
dislocations of Burgers vector b because 

(nb>* » n E b* 

Where = b, is the Burners vector of the elementary dislocation. 
Hence a giant dislocation of Burgers vector nb could burst into n 
component dislocations each with Burgers vector b, followed by an 
overall reduction in the energy of the system when the crack has 
passed through an active slip plane. The total number of 


dislocations in the pileup is given by the relationCSll : 

r a‘ 

n a = — (2.3.10) 

6 b 

Where, a is a correction factor for the case of dynamic 
dislocations and moving cracks. 

Using eqns (2.3.1) and (2.3.2), the strain energy density of the 
moving dislocation is given by the relation! 

u”* = — ^ — f P* + P* P* '! ^~(p P P P‘*'P P^"*" 4?^P* 

<4 2E v ^11 ^zz ^99^ E V it ^*2 ^aa ^ ttJ 2G ^tz 

(2.3.11) 

The plots of eqn(2.3.11) are shown in Fig. (2.9) for different 
normalized dislocation velocity. 
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2.3.4 STRAIN ENERGY OF AN INCLUSION* 

If a matrix volume is filled with a spherical inclusion 
of radiusCr*) and the body forces removed, then the matrix will 
relax, and the stress and strain fields produced are given by 
eqn(2.1.11>. The strain energy stored in the inclusion is given 
byC563* 

Ui « -4- <y .Av (2.3.12) 

it rr 

Where Av is the change in volume of the inclusion and equals to 
4fttfr**. Using eqn (2.1.11), eqn(2.3.12) becomes 






r* 


(2.3.13) 


2.3.5 INTERACTION ENERGY BETWEEN CRACK AND DISLOCATION* 


The interaction energy between a crack and dislocation 
is taken as the average of crack to dislocation energy (Ec d) 
and dislocation to crack interaction energy (Ed -♦ c) per unit 
length of crack fronts 

Ec — d » — ^ — C (2.3.14) 
The subscripts 'c' and 'd' refer to crack and dislocation, 
respectively. For volume element r d0 dr dz, interaction energy 
Ed c per unit length of crack front is given by following 
relation: 


Ed.c = ) r dr dG 

The is the dislocation stress components at crack tip 

i-j' 

is obtained by using the transpose of the rotation matrix 

and relation (2.1.15) and (2.2.6). The relations are : 

= cos*d + sin*G - sin2G 

11 11 zz tz 

s sin*G + cos*G o'! sin2G </^ 
zz 11 zz IZ 

cf z via * c I 

99 11 22 ^ . 


(2.3. 15) 
, which 
(2.1.14) 
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“ -j— sin2d(<y“^ — + cos2d <y 


(2.3.16) 


For a moving crack using eqns<2.3.2), (2.1.6) and (2.1.7), are 


given by following relations: 

c Kl z 


c Kl z 

.1° TT®i?[’.r *'('»* ^.a)] 

«' ’»)3 


c ^ Kl I-. , 

E Vanr 

c - 
^ a O 

22 

^ K l 

12 E TWr 12 


(2.3. 17) 


Similarly for stationary crack, using cqns(2.3.2), (2.1.9) and 

(2.1.10), c^. . are given by following relations : 

'tj 


c 


- u 

li E 

T 2 nr ^11 


c 

— [f (O) 

- V 

22 E 

rxnr k * * 


» 0 
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C 1 

+ V Ki 

(8) 

s = — 

12 

E Yznr ^tz 


^ ^ f <d) <2e3sl8> 

tz E rznr tz 

Th® cyf . and both are functions of two independent variables r 

t Jr tj 

and hence integral (2e3e 15 ) can be solved for moving and 

stationary crack separately with the help of Simpson's three point 
numerical integration formula for different normalized crack tip 
speed « 

Nowy the interaction energy Ec -► d per unit length of crack front 
is given by the relation : 

Ec d » (2.3*19) 

Where 'Wt' is the distance from crack tip to the boundary of medium 
For stationary cracky using eqn(2.1.17) and integrating (2.3.19) ^ 
it become: 

^ Ki b (Vr “ Vr) ^ •x 

Ec d « Tarr sxn8 cos 


For moving crack, using eqn(2.1.19) and integrating eqn(2.3.19), 


it becomes: 
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(2.3.21> 

The variation of interaction energy between crack and dislocation 
with & and for different normalized crack tip speeds have been 
shown in Fig. (2.10). 
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CHAPTER 3 

RESULTS AND DISCUSSIONS 

3.1 FORCE ON A DISLOCATION* 

Fig. (2.2) shows the variation of force F with B for 

static and moving cracks for different normalized crack tip speeds 

(A/c^>. The plots have two fold symmetry around the crack edge. 

For crack speeds less than about twenty percent of the shear wave 

speed c^ of the material, the force has a maKimum value along d = 
o 

70 . For crack speeds greater than 0a2c^^ however, this force has 

a maximum in a direction perpendicular to the direction of crack 

motion (d « 90^)* For crack speeds greater than or equal to 0.8c f 

a 

the curves have two maxima. The first maximum, which is of a 
greater magnitude, is at d * ‘?0°, for all speeds greater than 

0.8c^. The second smaller maximum has its position along B « 30° 

at crack speed of 0.8c^ and along B ~ 40° for all speeds in the 
range 0.90c^ < A < c^, c^ is Rayleigh wave velocity and it is 
approximately equal to 0.93c^. Rayleigh waves are characterized 
as waves which propagate parallel to the surface of a 
semi— inf inite solid and whose amplitude is greatest at the surface 
and decreases exponentially with distance from the surface. It is 
neither purely longitudinal nor purely transverse. At the 

Rayleigh wave velocity the force on the dislocation becomes 
infinite. In the range 0 < & < c^, the magnitude of the maximum 
increases with increasing crack speeds (see figures (2.2a) and 
(2.2b>). The second smaller maximum has its position along B = 
50° for all speeds in the range < A < c^, however, the 
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inaQnitude of the (naximum decreases with increasing crack speeds 
further (see fig.(2.2c>>. The force becomes attractive or 
repulsive depends on speeds of the crack tip and angular positions 
of the dislocation. The Table (2.1) shows the nature of the 
maximum force. 

Any emitted dislocation from crack tip or lattice 
dislocation moves away from crack tip under the repulsive force. 
However, a dislocation comes to rest at a point on the slip plane 
where the force on the dislocation is balanced by the lattice 
frictional force or by any other dislocation barrier such as 
inclusions. The dislocations pile up at this barrier. The 
piled-up dislocations exert image forces on the crack tip. The 
image force is always attractive and inversely proportional to the 
distance from the crack tipC321. This image force is responsible 
for the DFZ near a crack tipC21. The crack tip shielding provided 
by piled-up dislocations against applied stresses, reduces the 
local stress intensity factor at crack tipC303. The reduction in 
local stress intensity factor, may delay the crack propagat ion or 
the crack may deviate from its original path until piled-up 
dislocations cross the barrier. In the process the ductile to 
brittle transition point may shift towards the brittle side for 
the moment. If the Burgers vector of the emitted dislocation, had 
a component perpendicular to the crack front, then emitted 
dislocations blunt the crack tip. The local stress intensity 
factor also is reduced by blunting. Therefore, it is necessary to 
increase the shear stress to propagate the crack further. This is 
in analogy to the work hardening of materials. 
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In the range 0.9c < & < c the lattice dislocations 

Z K 

present at slip planes that intersect the crack front move towards 
the crack tip under the attractive force of crack at orientations 
such as d = 40°. With progressive crack growth, the movement of 
lattice dislocations is accompanied by plastic deformation near 
the crack tip region. The plastic deformation gives rise to 
blunting of crack tip and thus hinders the crack propagation. The 
movement of lattice dislocations at the crack tip during crack 
growth is therefore important in the case of ductile materials. 

3.2 EFFECT OF INCLUSION CW CRACK-DISLOCATION INTERACTEONt 
3.2.1 CRACK-INCLUSION INTERACTION: 

In Fig. <2. 3) and Fig. (2.4), the equipotenti al curves 
for static and moving cracks have been plotted. The two possible 
cases c > 0, and ^ < 0 are separately considered. 

Case I — If e: is positive, the interaction energy is lowest ahead 
of a crack tip, i.e. B = 0°, for static as well as moving cracks. 
In other words, the attractive force on the inclusion is highest 
ahead of a crack tip, i e d = 0°, for static as well as moving 
cracks. An inclusion with a greater elastic modulus such as many 
carbides in ot-Fe matrix produces a tangential compressive strain 
near the inclusion-matrix interface. Hence, under the radial 
attractive force, a crack could divert itself around the 
inclusion. For a static crack (Fig. (2.4)), the interaction energy 
increases from minimum at $ * 0 to maximum at B - ti for all 
equipotent ial curves. For a moving crack (Fig, (2.3) ) , the 
equipotential curves have two fold symmetry around the crack tip. 
The interaction energy increases from minimum at d = 0° to maxima 



36 


at 0 * n/2, and again decreases from maximum at $ = nfz to minimum 
at d = nr. For a static crack, the interaction energy is zero for 
d = n, however , for a moving crack the interaction energy never 
becomes zero for any angular position. 

Case II — If « is negative, the interaction energy is maximum 
ahead of the crack tip, t e O - 0°, for a static as well as for a 
moving crack. In other words, the repulsive force on the inclusion 
is highest for the position & » 0°. 

3.2.2 DISLOCATION-IHCLUSION INTERACnONt 

Case I - If « is positive, the interaction energy is negative for 
n < & < 2rt, hence an inclusion will attract a positive edge 
dislocation with its compression side lying on the side of the 
inclusion. The most favorable position for inclusion is at d * 
Zn/Zy where the interaction energy is minimum or attractive force 
is maximum. The inclusion will repel a positive edge dislocation 
with its tension side being on the inclusion for the positions in 
between O < 6 < n. 

Case II — If a? is negative, then from eqn(2.2.12>, it is clear 
that interaction energy is negative for, O < d < 2n, hence an 
inclusion causing a negative misfit strain will attract a positive 
edge dislocation with its tension side lying on the side of the 
inclusion. The most favorable position for inclusion is at 0 » 
n/Zj where the interaction energy is minimum or attractive force 
is maximum. The inclusion will repel a positive edge dislocation 
with its compression aide being on the inclusion for tz K & ^ 2tt. 

3.2.3 CRACK-DIS5LOCATION-INCLUSION SYSTEM* 

The quantitative effect of inclusion on 
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crack-dislocation intaraction is shown in Figures (2. 6^) • The 
Fig <2. 6a) shows the variation of force with B on the dislocation 
in the presence of the inclusion. The maximum force acting on the 
dislocation changes its position to B » 60° from B - 70° for a 

static crack. The magnitude of the force is also increased. For 
the static crack the curves are no more symmetric around the 
crack. For static crack the maximum in the range B - n to is 

at 280°. For speeds in the range 0.2c^ < k < 0.8c^, the curves 

are two fold symmetric around the crack and the magnitude of the 
force increases for all angular positions. Table (3.2) shows the 
angular position for the maximum force on the dislocation at 

different crack speeds and for ^ » 30°, 45°, 60°, 70°, taken 

around the inclusion. 

3.3 STRAIN ENERGY DENSITY OF CRACK: 

The minimum strain-energy-density cri terion (MSED) , 
proposed by SihC423 postulates that a crack propagates along a 
direction of minimum elastic strain energy density. Fig. (2.7) 
shows the variation of the strain energy density with B around 
a crack for different normalized crack speeds and at a fixed 
radial distance r from crack tip. For the static crack, the strain 
energy density increases from zero at G * tr to maximum at B ® 70 

and again to a minimum at B = 0°. For a moving crack, the strain 

energy density plots have two fold symmetry around the crack for 

all speeds. With increasing crack tip speeds, the magnitude of the 
strain energy density increases for all angular positions. Taking 
Sih's postulate into consideration it may be concluded from 
Fig.-(2.7) that the likely direction of crack propagation is along 
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the direction 8 = 0° in the velocity range 0 < i < 0.6c . At i = 

2 

0.8c^, curves have minimum strain energy density along the 
direction 8 — 60 , and hence a crack will propagate in this 
direction. At speeds greater than 0.8c^, the curves have minimum 
strain energy density along the direction 8 = 70°, and therefore 
the most likely direction for crack propagation. 

Fig, (2.8) shows the variation of total strain energy 
per unit length of crack front with normalized crack tip speeds. 
In the speeds 0 < & < c^(f ig. (2,8a) ) , the total strain energy 
increases continuously and variation is logarithmic in character 
because there is a net flux of energy into the crack tip. As 
speeds tend to the Rayleigh wave speed c^, the total strain energy 
of crack tends to infinity. This indicates that the crack becomes 
more and more unstable as crack tip speed increases and in the 
vicinity of Rayleigh wave speed, the crack is most unstable. The 
possible bifurcation of crack may take place to minimize the 
energy of the system. The increase instability of the crack with 
increasing velocity is due to the lower energy release rate from 
the crack tip. In the speed range c^ < i < c^, the total strain 
energy of crack decreases with increasing velocity from infinity 
as . shown in fig. (2.8b) because there is a net flux of energy out 
of the advancing crack tip and decreased instability of the crack. 
3.4 STRAIN ENERGY DENSITY OF A NOVING DISLOCATE ONt 

The Fig. (2.9) shows the variation of strain energy 
density with angle around a positive edge dislocation for 
different normalized dislocation velocities for a fixed radial 
distance r. The plots have two fold symmetry around the 
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dislocation. For ^ = o, the maximum of the strain energy density 

decreases with increasing velocity of the dislocation up to 0.4c . 

z 

At v/c^- the plot has maximum strain energy density at fi - 

40 and at v/c^* O.S, maximum located at = 60°. For fi - 90°, 

energy density increases with increasing velocity of 
the dislocation. 

3.S DISLOCATION 60IERATICa4 BY CRACKt 

Shear stresses from moving cracks can cause dislocation 
sources to operate on the slip planes that intersect the crack 
front. An inclusion present on the slip plane can also act as a 
source for dislocation loop formation etc., when an emitted 

dislocation encounters the inclusion and crosses it, leaving 
behind a dislocation loop around the inclusion. To obtain an 
equivalent number of dislocations to match the crack shear stress, 
we equate the strain energy density of a crack at an angular 
posi t ion ( r , ©> with the strain energy density of the moving 

dislocation at an angular positioned,/?). As an example, the 
number of dislocations have been calculated for dislocation 

velocities in the range O < v < 0.9c^ at normalized crack tip 

velocities 0.2, 0.6 and 0.8 and the results are plotted in 

Fig. (2.11) for & * 70°, (i = 0°, and r * d * 10/im. Fig. (2.11) 

shows that the number of edge dislocations equivalent to the crack 
shear stress increases with increasing dislocation velocity and 
variation is logarithmic in character for all crack velocities. 
However, with increasing crack velocity, the graphs move upwards. 
The number of dislocations increases with increasing crack 
velocity at all dislocation velocities. There is more pronounced 



40 


increase 

in the number 

of dislocations in 

the 

velocity range 

0 a 8c & 

< c^. The 

graph is restricted 

to 

the dislocation 

velocity 

less than the 

Rayleigh velocity c 

and 

the number of 


4 

dislocations is of the order 10 . 

Fig. (2.12) shows the distribution of the number of 
dislocations with normalized dislocation velocity and for 
different radial distance between crack and source for a fixed 
crack tip velocity(d » r » 5, 10, 20, 30, 50, 100pm, e * 

70°, ft * 0°, and k - 0.8c^) . The graphs moves downwards as the 
distance r increases. The number of dislocations decreases but 
the trend of graphs is the same for all radial distances. The 
Fig. (2.13) shows the distribution of number of dislocation with 
radial distance r and for Q * 70°, ft - 0°, d ® 10pm, 0.2, 

0.6, 0.8 and k - 0.8c . There is decrease in the number of 
dislocations with increase in distance. The graphs move upwards 
as dislocation velocity increases. The Fig. (2.14) shows the 

distribution of number dislocations with shear stress acting on 
the dislocation and for & * 70°, ft - 0°, d = 10pm, r = 10pm, and 

v/c^» 0.2, 0.6, 0.8. The plot shows the increase in number of 
dislocations with increase in shear stress as a function of the 
dislocation velocity. The Fig. (2.16) shows the variation of 
number of dislocations with ^ * 10pm, 

and V = k ~ 0.8c . The number of dislocations decreases with 
increasing the quantity o-^^r and variation is approximately 
hyperbol ic . 

3.6 DISLOCATIWf SUB-STRUCTURE* 

An acting dislocation - source producing n number of 
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dislocations at any tine, can have an equal number of positive and 
negative edge dislocations out of the total number n, with one 
half moving away from the crack tip under strong repulsive forces 
of the crack for B between to 85^ and for crack tip speeds in 
the range 0m2c^< i < 0«6c^(see Fig«2«2a>. These positive edge 
dislocations once formed relax into the matrix. The remaining n/z 
negative edge dislocations move towards the crack tip under strong 
attractive forces of the crack. This is reversed in the lower of 
the crack plane for B between 265°-355° and for crack tip speeds 
in the range 0.2c^< k < 0.6c^, the positive edge dislocations move 
towards the crack tip and negative edge dislocations move away 
from the crack tip and relax into the matrix. The movement of 
edge dislocations towards the crack tip results in many 
sub-structural changes ahead of crack tip and some of the 
interesting cases are considered in the ongoing discussion. 

Case I; Suppose there is no barrier for dislocation movement^ 
except the moving crack surface. The movement of dislocations is 
accompanied by plastic deformation. As the density of 
dislocations builds up at the crack tip, the plastic deformation 
is sufficient to slip or blunt the crack tip. The possible 
resulting crack geometry is shown in Fig. (2.17). The crack tip 
stress field promotes the generation of dislocations from a 
source, since it reduces the back stress on the source by pulling 
away the dislocations from the source. When blunting occurs, 
crack relaxes for a moment and then moves forwards. Crack tip 
blunting reduces the stress intensity factor and this results in 
slowing down the generation of dislocations at that particular 
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slip plane containing the dislocat^ion source. 

Case Ils Suppose there are barriers such as an inclusion in 
between the crack tip and the dislocation source, where the 
inclusion is harder than the matrix. The n/z dislocations pile up 
at the inclusion leading to the possible nucleation and growth of 
a new crack as shown in Fig. (2.18). Since the edge dislocations 
cannot climb like screw dislocations, the only option left for 
dislocations is to relax at the barrier by forming a new crack. 
The direction of the new crack is approximately at 90° to the slip 
plane of of dislocations. This case is possible in the case of 
relatively brittle materials. 

Case lilt If inclusion is harder and lies in a ductile matrix. 
The nucleation and growth of voids at the inclusion-matrix 
interface may be facilitated due to the build up of a tensile 
hydrostatic stress component from dislocations stresses in an 
already otherwise tensile field of the crack. The crack attracts 
the inclusion towards the crack tip and the attractive force on 
the dislocations also favors the displacement of inclusion towards 
the crack tip that results in building up of a 9ap between the 
inclusion and the matrix surfaces as shown in Fig. (2.19). This 
gap acts as a nucleation site for void growth. These mechanisms 
are some of the possibilities in the case of ductile fracture. 
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CHAPTER 4 

SUMMARY AND CONCLUSICWS 

1> The shear stress o*' associated with a moving crack has a 

12 ^ 

different spatial dependence than that of a static crack and 

has two fold symmetry around the crack tip. 

Si For crack speeds less than or equal to 0.2c an edge 

2 

dislocation experiences a maximum force on slip plane 
inclined, at d = 70°. For crack speeds in the range 0.2c^< A 

< 0.8c^, the maximum is along 9 - 90°. For crack speeds 

greater than or equal 0.9c^^ the curves have two maxima. The 

first maximum of greater magnitude is along & ~ 90° for all 

speeds greater than 0.8c^. The second smaller maximum is 

along O « 30° for crack speed 0.8c^, For crack speeds in the 

range 0.9c^ < i < c^ second maximum is located in the 

direction 9 = 40°. For crack speeds in the range 0 < A < c^, 

the magnitude of force increases with increasing crack speeds 

for all angular positions. For crack speeds in the range c^ 

o 

< A second maximum have position along 9 = 50 j but 

magnitude of force decreases with increasing crack speeds 
above c . 

R 

3:> The force on the dislocation tends to infinity at c^= 0.93c^, 
where c^ is the Rayleigh wave velocity. 

4> An inclusion with a positive misfit strain experiences a 
maximum attractive force along 0=0° from a moving as well 
as a static crack. If the elastic modulus of the inclusion 
is greater than that of the matrix, diverts the crack around 
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the inclusion. 

S> For a negative misfit strain, inclusion will attract a 
positive edge dislocation with its tension side lying on the 
side of the inclusion. 


For a positive misfit strain, inclusion will attract a 
positive edge dislocation with its compression side lying on 
the side of the inclusion. 


7> 

The presence 

of an inclusion changes 

the 

pos i t i on 

of 

the 


maximum force 

on the dislocation. The 

maKimum force 

on 

the 


dislocation changes its position to 9 = 

60° 

from & = 

0 

O 

N 

for 


static crack. 

With increasing value of 

the 

angle (^> 

around 


the inclusion 

the force curves becomes 

more 

f lat . 




For a static crack, strain energy density is zero at crack 

tip(d * n) and has a minimum in the direction & = 0°. For a 

moving crack, the strain energy density plots have two fold 

symmetry around the crack. The crack will propagate along 0 

- 0°, since the crack has minimum strain energy density along 

© * 0° for crack speeds in the range O < k < O.^c . At 0.8c^ 

z X 

strain energy density for crack is minimum along & ® 60° and 
for speeds greater than 0.8c^ curves have minimum strain 
energy density along & = 70°. 

9:> Total strain energy of crack increases with increasing crack 
speeds in the range 0 < k < c^ and tends ti infinity at 
Rayleigh speed c^. The instability of the crack increases 
with increasing crack speeds. At high speeds the crack 
bifurcates to minimize the energy of the system. At speeds 

greater than the c total strain energy of crack decreases 
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with increasing crack speeds, 

lOD The equivalent number of dislocations required to match the 

shear stress of a moving crack increases with increasing 

crack speeds, dislocation speeds and shear stress. The 

equivalent number of dislocation decreases with increasing 

radial distance from the crack tip. The number of 

dislocations required to match the shear stress is of the 
4 

order 10 and has a non-linear variation with the radial 
distance from the crack tip. 

115 The build up of dislocations at the crack tip results either 
in blunting of the crack tip when the crack surface acts as a 
barrier or leads to formation of a new crack in a direction 
perpendicular to the slip plane in the presence of an 
inclusion between the crack tip and the dislocation source 
or nucleation of a void at the inclusion and matrix interface 
depending upon the strength of the inclusion. 
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APPENDICES - A 

For small displacement a(t) <<< x and from Fig. (2.1) 
tan© = — X * r cos©, y = r sin© (A.l) 

C*)*''* (A. 2) 

Using eqn(2.1.3) and (A.l), eqn(A.2) becomes 
r* = jV cos© - a(t)]*+ <5* r* sin*© 

= r* cos*© + a*(t> - 2 r cos© a(t) + ©* r* sin*© (A. 3) 

Neglecting the term a (t) and 2 r cos© a(t) ,(A.3} becomes 
r* * r* £l - sin*© + <5* sin*© J 


or, 

r . = 

r ^1 - (1 - <5* > sin*©^*^* 

(A. 4) 

Now, 


©^ = tan (< ) 

(A. 5) 

Using 

eqn (2. 

1.3) , «qn (A. 5) becomes 




' X - a(t) 


X >>> 

a(t ) , 

tan©. * <5. — 

t X 


or, 


.tand , = <5 . t an© 

(A.©) 
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'5 ^ w 

! i.T., k; 


APPENDIX B 




aM22S1 


a. 


L*HOSPITAL*s THEOREMS Let ^(x) and ^(x> be the functions which are 
expansible by Taylor's Theorem in the neighborhood of x 
Also, ^<a) * O, and y/(s) - 0. Then 

- Ui» 


Lim 


'x-»a ^(x) ' x-+a yf' (x) 

In eqns(2.1.7) when velocity of crack tip & tend to zero all the 

0 


components of stress become the form of 


Hence this is an 


indeterminate form, so limit is obtained by using L 'Hospital's 
Theorem. Function D is common in all equations and this become 
zero, when 5. O. Differentiate eqn{2.1.8) w.r.t & and i 


becomes : 
dD 


d<5 


as 


= 4<5 


+ 4<5 


— - ^< 5 , (1 - 


as 


(B, 1) 


a& ^ ak ak 

Taking the limit of eqns<2.1.3) A ^ 0, we get 


Lim. rt <5 = Lim, , - 1 

k -> O t k ■¥ 0 z 

Differentiate the eqns<2.1,3) w.r.t A, we get 


ak 


as 


as 


k 


* X 


(B.2) 


{B.3) 


ak *"i ^ 

Substitute the values from equations (B.3) in eqn(B.l) and taking 
the limit as k tends to zero and simplified it becomes 


Lij 


at it 
ak 




r 1 ao ' 

1 = - 4 r ^ 

[A ak 

1 1 

1 

eqn (2.1.7) 

different iate 

as 

2<5 * 

fl + 2<5* - <5,) 

L * dA 

L 1 z-' 


V] 


(B.4) 


as 


ak 


2<5 
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- (1 . 4') (1 . 26* . 6*) 


2^ sin(^i/ 2 > 

2J '2 


r- 

as . ^ . 


2^ COS{dt/2> 


2 r ak 
1 


- ^(<5, 

ae 


+ 2S 6 
1 2 


26 6 

± 2 


COS (62/ 2) ^2 

2 y 
^ 2 


(B.5) 


Differentiate eqn<2«1.5) w.r.t. 4, uie get 




ae. 

tan© 

as. 

t 




ak 

1 + 6^, tan^€^ 

L 

ak 


Taking 

the 

limit of 

eqn (2.1*5) , 

we get 




6 , = 6 as 

t 

O 

t 



Differentiate eqn<2 

.1.9) w . r . t . 

4 and on arranging 

“ 3 / 2 

ae . 
1 

- = - 4 r 


f 4 sind 



ak 

^ L 

■=i i v 

^ c. J 

t 

J 

T ak ing 

the 

limit of 

eqn (2.1.9), 

we get 



(B,6) 


(B-7> 


<B.8) 


y.^lasSt-fO (B.9) 

^ t 

Put the value from eqns(B.3)^ (B.6> and (B-8) in eqn(B«5> and 
taking the limit, we get 

'-^^k of-r ] = ■ ^["V ■ ■V][^ ■ 3in36/2jcos6/2 

(B.IO) 

Divide eqn(B.lO) and Eqn<B.4), we get fii(6) for static crack; 


f 11(6) 




sind/2 sin3^/2l cos^/z 


(B.ll) 


Now from eqn(2.1.7) differentiate fzz component of stress w.r.t*4 
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df*2 


= - 46 


(1 * K) 


Zx cos iOt/z 

zJ ^ 


Cl - 


+ (1 + 6^y — 1- +4(6 

/; ^ 'i cos (02/2) _ r sin (02/2) *^2 

‘ tf4 ’''^2 * ' 2 yj-^ 

&y 

26 6 (02/2) ^2 

t°2 _ a / 2 „ 

2 >- 9h. 

* 2 


(B.12) 


Put the values from equat i ons (B . 3) , (B.6) and (B.8) in eqn(B.12) 


and taking the limit, uie get 


4+ ^]---[4r-4rl( 


1 + sind/2 sin3^/2jcos6?/ 2 


Divide the eqn(Bal3) and eqn<B«4), we get fzziO) for static crack 


f22(d) = |l + sin6^/2 sin3d/2j cos^/2 


CB. 14) 


Now from eqns(2.1.7) differentiate the fi2 component of stress 


w - r . t & 
&ftz 

a& 


COS(dl/2) 1 

2 


s i n ( 02 / 2 ) ^2 

2 r 6& 

* 2 


COS <02/ 2) ^2 


s i n ( 01 / 2 ) ^ ± 

2 r «k 

* 1 


] (‘ - <) * 


(01/ 2) sin (02/2) 


(B.15) 


Put the values from equat i ons (B . 6 > and (B.B) in eqn(B.15) and 


taking the limit we get 


r 1 6f 12 1 _ _ , r 1 _ 1 *1 

+ 0 [ i J " [ 2 2 J 


sin0/2 cos30/2 cos0/2 


(B. 16) 



Divide the eqn(B.16> and (B.4), we get f±zi&) component of 
for static crack: 


stress 


f 12(d) 


sin(d/2) cos (d/ 2 ) cos(3d/2> 


(B.17> 





1 » 1.12 '"0«5<S 0*0 0*5<b 1-12 


F cos© ^ 


Fi9.(2*2a) — The variation of force on a dislocation 

(F = O'' b) around a crack tip for a fixed radial 
tz 

distance in velocities range 0.0 < S. :< 0-8c^. 

o*' - Transformed shear stress on the glide plane of 
12 

d islocat ion 

b(Burger vector of dislocation) = 2.48A 
v(Poisson ratio) = 0.333 
GCShear modulus of matrix) = 73 Gpa 
KiOtress intensity factor) = 65 hpai/m 
r(Radial distance) = lO/jm 




Fig. (2.2b> — The variation of force on a dislocation 

(F = O'' b ) around a crack tip for a fixed radial 

distance in velocities range 0.90c < k < 0.92c . 

2 2 

o' = Transformed shear stress on the glide plane of 

12 ^ K 

dislocation 

b = 2.48A, V = 0.333^ G = 73 Gpa, Ki = 65 MpaVm, 
r = 10/L^m 






F cos6^ > 


Fig, (2,2c) — The variation of force on a dislocation 
(F = b) around a crack tip for a fixed radial 
distance in velocities range 0,94c^ < k < c^. 

O'' = Transformed shear stress on the glide plane of 

dislocat ion 

b = 2. 48 A, V = 0.333, G - 73 Gpa, Ki = 65 lipaVm, 


r = iO^^m 


giiiniiijtminiTjminTrr[nnimr|iii^'YfC”Ppfi''Titijmiimi|TirminjtrtniTrr| 







Fig. (2.3b) - Equipatent i al contours of constant value of 
D ^ <5^)cos (^ 1 / 2 ) at crack tip velocity 

&, = 0,4c^ for moving crack and inclusion interaction* 



Fig.<2-3c) - Equipotent i al contours of constant value of 
Q - p” (l ^ 2 ) crack tip velocity 


§L = 0.6c for moving crack and inclusion interaction, 
2 ^ 




Fig. (2.4) - Equipotent i al contours of constant value of 
1 , 

“yp cos C 0 / 2 ) for a static crack and inclusion 


interaction . 
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- 1.8 - 1.2 - 0.6 0*0 0.6 1.2 l-B 


F cos^* > 

Fig. (2.6a) - The variation of force on a dislocation 

T 

(F = O' b) around a crack tip for a fixed radial 
±2 

distance in velocities range 0.0 < ^ < 0.8c . 

a = 0*^0 
12 12 12 
T 

a = Total shear stress on the dislocation 
1 2 

o^ = Transformed shear stress on the alide plane of 

±2 - r ' 

dislocation 

o^^= Transformed inclusion shear stress on glide plane 
of dislocation for = 30^ and r‘’/r'= O.B 
b = 2.48A, V - 0.333, G = 73 Gpa, Ki = 65 MpaVni, 

/j = 156 Gpa, r = lO^m, r® = 4.2A. 



o / 



Fig, (2.6b) - The variation of force on a dislocation 

(F = a'*' b ) around a crack tip for a fixed radial 
12 

distance in velocities range 0,0 < k < 0,8c , 

2 

Y I 

ty = <5* O* 

12 12 12 

!= Transformed inclusion shear stress on glide plane 
12 ^ 

of dislocation for (p = 45 and po/r'^ 0.8 
b = 2. 48 A, V = 0.333, G = 73 Gpa, Ki = 65 MpaVm, 
p = 156 Gpa, r = 10/jm, r® = 4,2A, 



. a& 



- l.B - 1,2 -" 0.6 0.0 0.6 1-2 1-8 

F CQS^ > 


Fig, (2-6c) The variation of force on a dislocation 

T 

(F = around a crack tip for a fixed radial 

distance in velocities range 0-0 < i < 0.8c 

T , I ^ 

<y = cy' 4- cy 

j 12 12 12 

^ 12 ” Transformed inclusion shear stress on glide plane 
of dislocation for <p - 60° and 0.8 

b = 2,48A, 'U = 0.333^ G = 73 Gpa, Ki = 65 MpaVm, 
fj = 156 Gpa, r = lO^m, = 


4,2&, 



sin^ 



FiQ-{2.6d) ~ The variation of force on a dislocation 

(F = O' b) around a crack tip for a fixed radial 
12 

distance in velocities range 0.0 < k < 0.8c . 

2 

= Transformed inclusion shear stress on glide plans 
i2 O 

of dislocation for <p = 70 and r^/r'= 0,8 
b = 2. 48 A, V = 0.333, G = 73 Gpa, Ki = 65 MpaVm, 

^ = 156 Gpa, r = lO^m, r® ~ 4,2A. 
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Fig. (2-7) ~ The 

cos^xlO^ ► 

distribution of 

the elastic 

strain 

energy density at 

the crack tip for 

values of 

B 

varying 

between 0^ and 

360^ and for normalized 

crack tip 

velocities 0-0- 0 

.2, 0.4, 0.6, 0.8. 
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sin© 



-1.2 -0.8 -0.4 0*0 0.4 0.8 1.2 


u”? cos© 1 

a 

Fiq.( 2,9) - The distribution of the elastic strain 
energy density around the dislocation for value of d 
varying between and 3^0^ and for normalized 

dislocation velocities 0.2, 0.4, 0.6, 0.8. 



Fig. (2.10) — The variation of interaction energy E 

c-d 

between crack and dislocation with d for a fixed radial 
distance and for normalized crack tip velocities 0.2, 
0.4, 0.6, 0.8. 















of D j s 1 ocat i on (n )x 10e4 
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Fi9(2.16) - Showing the distribution of dislocations 

with (O'' r) for & = v = 0.8c . 

12 2 





Fi 9 «( 2 . 19 ) - Th© formation of void at inclusion and 
ductile matrix interface^ when inclusion is harder than 
the matrix. The hydrostatic stress field of the piled up 
is predominent ly tensile in the upper half of the glide 
plane and lies in the tensile field of crack. 


TABLE NO . 2 • 1 “ The nature of maximum force acting on dislocation 
at different angular positions and normalized crack 

tip speeds. 








TABLE NO-2.2 - The nature of force at different angular position 
for crack-dislocat ion- inc lus ion system . 

(+) - Repusive Force 
(-) - Attractive force 











